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We calculate how the relic density of dark matter particles is altered when their annihilation is 
enhanced by the Sommerfeld mechanism due to a Yukawa interaction between the annihilating par- 
ticles. Maintaining a dark matter abundance consistent with current observational bounds requires 
the normalization of the s-wave annihilation cross section to be decreased compared to a model with- 
out enhancement. The level of suppression depends on the specific parameters of the particle model, 
with the kinetic decoupling temperature having the most effect. We find that the cross section can 
be reduced by as much as an order of magnitude for extreme cases. We also compute the fi— type 
distortion of the CMB energy spectrum caused by energy injection from such Sommerfeld-enhanced 
annihilation. Our results indicate that in the vicinity of resonances, associated with bound states, 
distortions can be large enough to be excluded by the upper limit \fj,\ < 9.0 x 10~ 5 found by the 
COBE/FIRAS experiment. 
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I. INTRODUCTION 

If dark matter annihilates, the byproducts of the 
annihilation (positrons, neutrinos, gamma-rays, etc.) 
can leave non-gravitational signatures that, if observed, 
would be crucial for clarifying the nature of dark matter. 

In recent years, a number of observations have high- 
lighted anomalies that might be explained by invoking 
dark matter annihilation in our Galactic halo. Among 
them are: (i) an anomalous abundance of positrons in 
cosmic rays above 10 GeV according to the PAMELA 
experiment [l[ confirming and extending previous mea- 
surements by experiments such as AMS-01 2]; (ii) an 
excess of microwave emission from the galactic center 
as measured by the WMAP experiment, known as the 
"WMAP haze" 3 (iii) the apparent excess of diffuse 
galactic gamma-rays with energies above lGeV inferred 
from observations by the EGRET satellite Q ; (iv) anal- 
yses on the balloon experiments ATIC and PPB-BETS 
p, ||| have reported an excess in the total flux of electrons 
and positrons in cosmic rays. For the last two anoma- 
lies, we note that recent observations by FERMI seem 
inconsistent with the claimed gamma-ray excess in the 
EGRET data 0] , and that the excess in the electron and 
positron flux is smaller than previously thought Q . The 
latter is actually consistent with a modified cosmic ray 
propagation model that does not require additional pri- 
mary sources of electrons and positrons [§] . 

Although other astrophysical sources could explain 
these anomalies (see for example for an expla- 

nation for the PAMELA excess based on particle accel- 
eration by pulsars, and [T3, [l4| for one based on super- 
nova remnants), dark matter annihilation offers an at- 
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tractive solution. Large annihilation rates are needed, 
however, to explain the observations. In particular for 
the PAMELA data, the required annihilation rate is typ- 
ically a few orders of magnitude larger than the value 
obtained by assuming the standard cross section inferred 
from the observed abundance of dark matter together 
with a smooth local distribution of dark matter [Tlj . 
Thus, an additional hypothesis is needed to boost the 
annihilation rate to the required levels and this must not 
change the present-day abundance of dark matter. 

Such a boost is difficult to obtain from the effects of 
substructures in the local dark matter distribution. Re- 
cent numerical simulations predict this to be remarkably 
smooth [l6j]. A detailed analysis of the im pac t of sub- 
structure on the production of positrons by 17] came to 
a similar conclusion. The possibility of a nearby "spike" 
of dark matter produced by an intermediate mass black 
hole seems also a priori implausible [l8j . 

An alternative that has produced a plethora of papers 
in recent years is that of a Sommerfeld enhancement to 
the cross section produced by the mutual interaction of 
the annihilating dark matter particles. Such an interac- 
tion could be produced by a force carrier which might be 
any of the standard model weak force gauge bosons [l9[ 
or a new force carrier [2(j. For certain values of the pa- 
rameters of these models, the enhancement is easily large 
enough to boost the cross section to the required values. 

However, a large cross section has a significant im- 
pact on other observables and may violate other con- 
straints. For instance, 2l| showed that for the case where 
the cross section increases as 1/v (a particular case of 
Sommerfeld-enhancement models), where v is the rela- 
tive velocity of the annihilating particles, there are se- 
vere constraints from measurements of the diffuse extra- 
galactic gamma-ray background radiation and from CMB 
constraints on ionization and heating of the intergalactic 
medium (IGM) by annihilation in the first generation of 
halos. In this scenario, the boost factors required to fit 
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the above anomalies would be inconsistent with current 
constraints. This problem is avoided, however, by more 
general cases of the Sommerfeld enhancement where the 
effect saturates at low velocities [2(| ■ 

Recently, [22| and [23| analyzed constraints on the an- 
nihilation cross section from perturbations to the CMB 
angular power spectra resulting from heating and ion- 
ization of the photon-baryon plasma at recombination. 
They found interesting upper limits to models with Som- 
merfeld enhancement (see fig. 5 of (22j) that already rule 
out some extreme cases. 

At higher redshifts, the effects of the Sommerfeld en- 
hancement have scarcely been treated. For example, the 
annihilation of dark matter impacts the predictions from 
big-bang nucleosynthesis on the abundance of light ele- 
ments (e.g. [13] )■ This could put constraints on certain 
models with Sommerfeld enhancement. However, this 
has only been studied in passing (2|| |26| . 

Also, the abundance of dark matter today is commonly 
assumed to be unaltered by this effect, apparently be- 
cause the typical dark matter particle velocities at freeze- 
out are very large ma king the enhancement very close to 
one at that epoch [2(J[27[- As we will show, this reason- 
ing is flawed. 

The thermodynamic equilibrium between matter and 
radiation in the early Universe would be perturbed by 
energy released during a certain process, dark mat- 
ter annihilation for example. This equilibrium tends 
to be restored by different interaction mechanisms: 
Compton scattering, double Compton emission and 
bremsstrahlung radiation. The efficiency of these to fully 
restore equilibrium varies with redshift. For z > 2 x f 6 
they are efficient enough to restore distortions in the en- 
ergy spectrum and thus, the photon distribution is that 
of a black body with a slightly higher temperature than 
the one in the case of no energy injection. For lower red- 
shifts these mechanisms can not restore the black body 
spectrum. In particular, for 5.1 x 10 4 <z<2x 10 6 , the 
spectrum is perturbed into a Bose-Einstein distribution 
with a chemical potential fj, [281 ] . 

In this paper, we revisit the impact of the Sommerfeld 
enhancement on the relic particle abundance and we show 
that its effect is not negligible. We also study, for the first 
time, fj,— type distortions of the CMB spectrum due to 
energy deposition by dark matter annihilation in models 
with Sommerfeld enhancement. 

The paper is organized as follows. In section 2 we 
summarize the Sommerfeld enhancement and describe 
how we include it in our calculations. The relic density 
calculation is set out in detail in section 3. In section 
4, the [i— type distortion to the CMB from annihilation 
with Sommerfeld enhancement is calculated. Finally we 
present a summary and our conclusions in section 5. 



II. THE SOMMERFELD ENHANCEMENT 

The Sommerfeld enhancement of dark matter anni- 
hilation is a nonrelativistic quantum effect occurring 
when annihilating particles interact through a potential 
of some sort. If their kinetic energy is low enough, their 
wave function is distorted, and a significant enhancement 
(or suppression if the force is re pul sive) to the annihila- 
tion cross section occurs [ljl [2(Jl2!|. The product of the 
relative velocity times the annihilation cross section will 
therefore be boosted: av = 5(<7v)o, where (ctw)o is the 
standard product of the tree level cross section times the 
relative velocity and S is the so called Sommerfeld boost. 

The simplest interaction is that for an attractive 
Yukawa potential with coupling strength a c mediated by 
a scalar boson of mass (more complicated models ex- 
ist where there are multiple force carriers (20[). In this 
case, the Sommerfeld enhancement can be computed by 
solving the radial Schrodinger equation for s-wave anni- 
hilation in the non-relativistic limit (we use natural units 
ft = c= 1): 

+ ^e'^(r) - -fhcfm (1) 

where m x and f3 are the mass and velocity of the dark 
matter particle, respectively, and \l/(r) is the reduced 
two-body wave function. Eq. ([T]) should be solved with 
the boundary condition: d^/dr = im x f3^ as r — > oo. 
The Sommerfeld boost S is simply given by: 



In the limit of a massless force carrier (to^ — > 0), the 
potential reduces to a Coulomb interaction and Eq. (p} 
can be solved analytically yielding a value of S that is 
independent of the dark matter particle mass: 

S=^(l-e-™ ^ (3) 

In this limit, and for small velocities (/3 <C 7ra c ): S ~ 
Tra c /f3. This is the reason why the Sommerfeld enhance- 
ment is usually associated with a "l/u" enhancement. 

In the more general case, where ^ 0, the 
behavior is no longer valid for very small velocities be- 
cause the Yukawa interaction has a finite range. This 
can also be shown by expanding the exponential term in 
Eq. ([T]) in powers of m^r, then the necessary condition 
to recover the Coulomb-like interaction solution (Eq. [3]) 
is p 2 3> ot c m,p/m x fl9| , which, for fixed values of a c and 
m^/m x , will not be fulfilled for arbitrarily low veloci- 
ties. In the opposite regime, for 1 -C Q c m^/m x , a se- 
ries of resonances appear for specific values of a c m x /m^. 
These resonances are associated with bound states, the 
so-called WIMPonium [3(1 |3l|. In them, the enhance- 
ment greatly increases with smaller relative velocity, as 
S ~ I//? 2 . For even lower velocities and for fixed values 
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of a c and m^/m x , the enhancement saturates reaching 
a value that can be estimated by solving Eq. ([T]) with 
(3 = 0. Away from the resonances, an order of magnitude 
estimate of this value is S max ~ 6a c m x /rn j,, reached for 
a threshold velocity of (3 ~ 0.5m r f>/m x [19j| . 

On the other hand, for large velocities (j3 ^> 7ra c ), the 
Coulomb approximation is valid and the enhancement is 
effectively reduced to values very close to 1. 

To include in detail the influence of the Sommerfeld 
enhancement in our study we need to to solve numerically 
the Schrodinger equation. 



III. RELIC DENSITY CALCULATION 

The evolution of the phase space distribution per parti- 
cle of a given species is given by the Boltzmann equation. 
For Majorana fermions that are stable, for example, neu- 
tralinos if they are the Lightest Supersymmetric Particle, 
the Boltzmann equation is [e.g. |32j: 



dt 



(av) (n x - (n^f) ( 4 ) 



where H = a/a is the Hubble parameter, (av) is the 
thermally averaged product of the total pair annihilation 
cross section times the M0ller velocity and is the 
equilibrium dark matter number density. In the remain- 
der of the paper we assume that dark matter is made of 
neutralinos. 

We solve Eq. (O by using the standard freeze-out ap- 
proximation [e.g. [32J . The freeze-out epoch is defined as 
the time where the annihilation rate is equal to the ex- 
pansion rate of the Universe (r ~ H). In the standard 
case where (av) is velocity independent, the annihilation 
process becomes subdominant after freeze-out (chemical 
decoupling) and the "relic" abundance of particles is 
roughly fixed [45j. 

We define A = X — Xeq, where X and Xeq are the 
comoving number density out of and at equilibrium, then 
we can write Eq. (U|) as: 



dA 
~d7 



dX EQ 
dt 



(av) A (2X EQ + A) (5) 



Since before freeze-out the number density follows closely 
the equilibrium solution then: dA/ dt <C dX E Q / dt, and 
at freeze-out: A(tf) m cX E Q(tf), where c is a constant of 
order unity 46] . Under these approximations, the value of 
the freeze-out time can be found by solving the following 
equation: 



d\nX E Q 



dt 



-3 



(av) f c(2 + c)X EQ (t f ) (6) 



It is convenient to change the variable X to the ratio 
Y = n x /s of the number density to the total entropy 
density of the Universe. The entropy is dominated by the 
contribution of relativistic particles s = (2n 2 /45)g,sT 3 , 
where g*s are the effective degrees of freedom for the 



entropy density, and T is the radiation temperature. In 
thermal equilibrium, the entropy per comoving volume is 
conserved sa 3 = s ~ 2918 cnT 3 (taking T = 2.728°K 
and g*s(T ) = 3.915). 

It is also useful to change the time variable to x = 
m/T using the equation: t = 0.301gt 1/,2 mp t (x/m) 2 = 
G,(m,T)x 2 /2, where m P i = 1.221 x 10 19 GeV is the 
Planck mass and <?* are the effective degrees of freedom 
for the total energy density of the Universe, which are 
generally a function of temperature. We have taken into 
account the temperature dependence of in our calcula- 
tions, whose values vary from ~ 100 at T ~ 100 GeV 
to ~ 3 at T ~ 10 keV. This equation is valid during 
the radiation-dominated era. If we neglect the variations 
in g*(T) around the freeze-out temperature [47|| . we can 
finally write Eq. ^ as: 



dlnFj 



EQ 



dt 



- ~G*{m, T f )s f (av) f c(2 + c)Y EQ (xf)x f 

In the non-relativistic limit, the equilibrium solution 
takes a simple form: Yeq(x) = 0.145(g/g*s)x 3 ^ 2 e~ x , where 
g is the degeneracy factor for the particle species (g = 2 
for neutralinos). Therefore, xt is simply given by the 
solution to the implicit equation: 



e Xi - 0.038 



<J 



,3* /2 (T/), 



mpi c(2 + c) 




av) f (8) 



where we have used the approximation 3x * 1 /2 <C 1 , since 
we expect Xf > 20. 

In the non-relativistic limit, the thermal average (av) 
reduces to an average over a Maxwell-Boltzmann distri- 
bution function I32f: 



(av) 



2tt 1 /2 



(av) /3 2 e-^ 2 / 4 d/3 



(9) 



where v is the relative velocity of the annihilating parti- 
cles. The expansion of (av) in powers of j3 leads to the 
usual expansion of the non-relativistic thermal average 
in powers of x, where the first term of the expansion cor- 
responds to the constant s-wave annihilation term (av) s . 
For the case when the cross section is enhanced by the 
Sommerfeld mechanism the s-wave annihilation thermal 
average is given by: 

M = Ms {^[ S ^ /3 2 ^ 2/ V) = {™) s S(x) 

where S(x) is the thermally averaged Sommerfeld en- 
hancement. The thermal average maintains the quali- 
tative behavior of S with velocity. For example, when 
S(j3) oc S(x) oc x 1 / 2 oc a~ l , where a v is the ve- 

locity dispersion of the dark matter particles [Hj]. We 
have omitted in the notation of S(x) the dependence 
on the parameters of the Yukawa interaction; in general 
S = S(x,a c , m^jm^). 
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At freeze-out, although the dark matter particles are 
already non-relativistic, the typical velocities are still 
very large, and we can safely assume S(j3) ~ 1, i.e., 
(av) f ~ (crv) s . We have checked that this is indeed 
a good approximation by solving fully the Schrodinger 
equation, there is no significant change in the results that 
follow. 

The dark matter relic density can be calculated by 
solving the Boltzmann equation in the late times regime 
(i > tj) where the number density at equilibrium n% Q is 
considerably lower than the out-of equilibrium solution 
and can therefore be dropped from the Boltzmann equa- 
tion. Thus, Eq. @ can be written in terms of the ratio 
Y = n x /s and the variable x = m x /T as [e.g. l32l|: 



x 1/2 

9 



(x)S(x) 



Y( Xo ) Y( Xf ) • V 45 mxmpi {<Xv) sJ Xf x 2 dx 

(11) 

where we have used Eq. (fTUf) , and the limits of integration 
corresponding to freeze-out Xf = m x /Tf and present-day 
values xq — m x /To. The value of Y(xq) is connected to 
the value of the ratio of the dark matter density to the 
critical density today (il x ,o = m x n x,o/ 

Per it, ) ■ 



8vr 



ZM 2 pl {H Q /h) 2 
~ 2.757 x 10 8 



m x s{x )Y(x ) 



GeV 



Y(x ) 



(12) 



The integral in the r.h.s of Eq. (fTTj) can be divided in 
two. Before kinetic decoupling, the dark matter particles 
are still coupled to the photon-baryon plasma through 
scattering with standard model particles and therefore 
we can treat the temperatures of the radiation and of 
the neutralinos interchangeably (T x = T) for x < Xkd — 
m x /Tkd, where Tkd is the temperature at kinetic decou- 
pling. 

After kinetic decoupling, the temperature of the dark 
matter particles drops as a~ 2 instead of a -1 , as the 
temperature of radiation does. Thus, for this second 
part of the integral in Eq. (fTTj) . one must be careful in 
the definition of the variable x in the thermal average 
(Eq. []JJ and in the Boltzmann equation (Eq. [TT]) . The 
thermal average uses a Maxwell-Boltzmann distribution 
function with the temperature of the dark matter "gas" , 
whereas Eq. (111[) was written as a function of x — m x /T, 
where T is the radiation temperature. We can easily ac- 
count for this by evaluating S(x x ) in Eq. (fTUf as a func- 
tion of the radiation temperature by using the relation: 
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FIG. 1: Values of the normalization of the thermally aver- 
aged cross section (crv) s that give the proper relic density 
h 2 = 0.1143 in the parameter space (m r f,/m x , a c ) for 



100 GeV and T kd 



MeV. The values are given 



in a color scale shown on the right hand side in units of 
1 x 10" 26 cm 3 s _1 . 



before matter-radiation equality. Consider, for example, 
a 100 GeV neutralino. Since z^~4x 10 3 , and because 
T x ~ {zE/zkd) 2 with Zkd ~ 4 x 10 10 (the exact value of 
Zkd depends on the specific supersymmetric model, this 
is a typical value for a neutralino of this mass), we have 
a dark matter velocity dispersion of the order of 10~ 8 
by z = ze, which is sufficiently low to reach saturation, 
unless the combination of a c and /m x is fine-tuned to 
be close to a resonance. 

In summary, a proper calculation of the relic density 
needs to include the effects of kinetic decoupling. We 
do so by choosing a model with neutralino dark matter, 
and by using the results of (35|, who made a scan of 
the parameter space in mSUGRA (minimal supergravity) 
and more general MSSM models, and give a range for the 
values of the kinetic decoupling temperature as a function 
of neutralino mass (see fig. 2 of their paper). Typical 
values for m x € (100, 5 x 10 3 ) GeV are in the range Xkd S 
(200,2 x 10 4 ^ 



x 2 /x kd . As pointed out by [3J], in this regime For a fixed value of c^/^ we can solve Eqf ^ jgj and 



1/2 



oc x. 



the 



(x > Xkd), for the case where S(x x ) oc x 
relic density decays logarithmically: Y oc 1/lnx. In sec- 
tion 2 we mentioned that for particular values of a c and 
m ( j ! /m x , the Sommerfeld enhancement has resonances 
that for very low velocities produce very large values of S, 



in these cases, S(x x ) oc x x oc 
decreases as 1/x. 

However, these large suppressions of the relic density 
are eventually cut-off by the saturation of S, which occurs 



pip simultaneously to get the values of Xt and (<tv) s 
that give the correct relic abundance. The value of 
the present-day dark matter density according to the 
WMAP-5 year data, at the 2a level, is fl x h 2 = 0.1143 ± 
0.0068 [36j. Fig. [JJ shows the value of (av) s consistent 
thus, the relic density with Q, x h 2 — 0.1143 for a scan in the parameter space 



100 GeV and T, 



kd 



8 MeV. 



(m^/m x , a c ) with 
The values of the thermally averaged cross section are 
color-coded in units of 1 x 10 _26 cm 3 s~ 1 . For compari- 
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<0~V >T M =150 10~ 3 Gev/ <<™ >T m =4 10" 3 GeV 




FIG. 2: Map ratio of the value of the cross section giving 
the correct relic density Q x h 2 = 0.1143 with different kinetic 
decoupling temperatures, T^d = 150 MeV and Tkd = 4 MeV 
for a 100 GeV neutralino. 



son, in the absence of any Sommerfeld enhancement the 
required cross section is ~ 2.4 x 10~ 26 cm 3 s _1 . 

Since the typical velocities of the dark matter parti- 
cles between Xf and Xkd are still very high, the Coulomb 
regime is valid in this range for most of the parameter 
space (recall that the Yukawa part of the potential can 
be ignored when 1 ^> a c m^/m x ) 7 thus, if we cut the 
integral in the r.h.s of Eq. (fTTj) at x = Xkd> the relic den- 
sity calculation would depend only in the strength of the 
coupling a c . This dependence can be seen in the overall 
trend from left to right in Fig. Q] However, after kinetic 
decoupling, the WIMPs get colder much faster than the 
radiation and the Sommerfeld enhancement depends not 
only on a c but also on m^/m x . This is particularly sig- 
nificant near resonances where the enhancement is large 
enough to suppress the dark matter abundance by a fac- 
tor of a few, therefore, the cross section normalization 
needs to be suppressed by this same factor in order to be 
consistent with the observed relic density. 

As shown by the discussion in the previous paragraph, 
the particular value of the kinetic decoupling tempera- 
ture has an impact on Fig.[TJ The closer this temperature 
is to the freeze-out temperature, the stronger the impact 
will be. To appreciate this better, we show in Fig. [2] the 
ratio of two maps analogous to the map in Fig. [T]but for 
two extreme values of the kinetic decoupling temperature 
for m x = 100 GeV: T kd = 150 MeV and T kd = 4 MeV. 
For large values of the coupling strength and low values of 
m ( j 1 /m Xl and particularly near resonances, the cross sec- 
tion at freeze-out needs to be lower for a higher kinetic 
decoupling temperature to produce the correct relic den- 



sity. From the figure we see that uncertainties in the ki- 
netic decoupling temperature introduce changes to {av)s 
at the most of order 30% for a 100 GeV neutralino. 

Since current uncertainties on the abundance of dark 
matter are in the percent level, the impact of the partic- 
ular value of Vt x h 2 on Fig. [T]is not so large. By exploring 
changes in fl x h 2 within the 2a bounds we found changes 
in the values of (av)s of 20% at the most. 

These results are also not very sensitive to changes in 
the neutralino mass, because the relic density is nearly 
independent of it. This can be seen by looking at Eqs.tfTlT- 
IT21) . If we neglect the term 1/Y(xt), then, for a fixed 
value of (<jv)s, £l x h 2 depends on the neutralino mass just 
through the integral in the r.h.s of Eq. (ITT1) . Moreover, 
since the Sommerfeld enhancement only depends on the 
neutralino mass through the ratio m^/m x , a change on 
m x can always be compensated by choosing a different 
mass for the force carrier. The remaining dependence is 
through the values of Xf and Xkd ■ 

Of course, to get the precise impact on the value of the 
cross section required for a particular particle model, it 
is necessary to do the calculation for that model. Un- 
certainties on the kinetic decoupling temperature and on 
the abundance of dark matter will have a minor-to-mild 
impact on the quantitative results. However, the general 
behavior and order of magnitude of the expected values 
of {av)s can be directly read off from Fig. Q] 

In summary, the relevant difference between the relic 
density calculation for the standard case where (crv) is 
constant and the case we study here is given by the on- 
set of the Sommerfeld enhancement after kinetic decou- 
pling. Before this phase, the velocities of dark matter 
particles are too large for the enhancement to be rele- 
vant, but once the decoupling happens the typical veloc- 
ities decrease quickly since the temperature of the dark 
matter fluid drops as a~ 2 . This causes a rapid increase 
on the enhancement values and the annihilation process 
takes a relevant role once more decreasing the dark mat- 
ter density. On the contrary, the comoving density re- 
mains constant after freeze-out for the standard case. 
The aforementioned effect continues to be important un- 
til the enhancement saturates fixing the value of (av) 
to a maximum. Afterwards, the annihilation rate be- 
comes subdominant compared to the expansion rate of 
the Universe and thus, the comoving dark matter den- 
sity remains fixed. 

Owing to this difference, the value of the annihila- 
tion cross section before the freeze-out epoch needs to be 
smaller for the case with Sommerfeld enhancement than 
for the standard one in order to produce the observed 
dark matter density. With a lower annihilation rate, the 
dark matter particles go through chemical decoupling at 
earlier times, and therefore the dark matter density is 
higher at these epoch than in the standard case. This 
higher value compensates the subsequent decrease of the 
density due to the Sommerfeld enhancement. In this way, 
the final relic density can be consistent with the expected 
value for the abundance of dark matter today. 
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IV. /i-TYPE DISTORTION OF THE CMB 
SPECTRUM 



Measurements of the CMB spectrum by the FIRAS in- 
strument onboard COBE have shown that is nearly that 
of a perfect blackbody with temperature Tq ~ 2.728° K. 
Deviations from a pure black-body spectrum are, how- 
ever, expected if there is any energy input in the early 
Universe. If this energy is injected between 5.1 x 10 4 < 
z < 2 x 10 6 (12 eV< T <470 eV), the spectral distortion 
is a Bose-Einstein \x— type distortion where /i is the chem- 
ical potential [28[ . If the number of photons injected dur- 
ing the energy release is small, compared to the number 
of photons in the radiation plasma, then /i ~ 1.4<5/3 7 //9 7 , 
where <5p 7 is the injected energy and /9 7 is the energy 
density of the CMB photons. 

The current observational limit for the //-type distor- 
tion at the 95% confidence level is \fi\ < 9.0 x 1CT 5 [jtfj]- 
In principle, this limit could already have been improved 
by nearly two orders of magnitude by now given recent 
technological advances [38}. 

In the case of energy deposited by annihilation of neu- 
tralinos, the p— type distortion was studied by [39| for 
the case of s-wave and p-wave cross sections. More re- 
cently it was also mentioned by [io] |. The value of /x is 
given by: 



/' : 



lA** = 1.4 

Pi 



El 

p-, 



dt = 1.4 



fm x (av) n\ 



(13) 

where, p 7i o is the present-day energy density of the CMB 
(fi 7 ft, 2 ~ 2.47 x 10~ 5 ), t\ — ti is the time interval corre- 
sponding to the energy injection and / is the the effi- 
ciency in which the injected energy is transformed into 
heat [48| . In principle, / depends on the channel of anni- 
hilation and on time. However at the relevant redshifts 
(z *3> 1000), / is basically given by the annihilation chan- 
nel [23j]. For electrons and photons / ~ 1, but for anni- 
hilation into r's for example, some of the energy is lost 
in neutrinos. Nevertheless, for all relevant channels, [23| 
found that / > 0.25 for z > 2500. We will adopt / = 1, 
our results can be reinterpreted easily for a different value 
of/. 

Since the relevant redshift range of energy injection is 
in the radiation dominated era and after kinetic decou- 
pling, we can write t = t kd {l+z kd ) 2 /(1+z) 2 = t kd x x /x x d . 
It is also convenient to write Eq. (|T3l in terms of the vari- 
able x x — m x /T x : 



P = 1-4/ 



tkd 



where: 



m xPifiJ \ x kd{^ + z k d) 4 



S(x x )x x 1 p 2 x {x x )dx x 



m x T kd 



„1,2 



(14) 



(15) 



FIG. 3: Scan of the parameter space (m$/m x , a c ) with the 
expected values of the p-type distortion to the CMB spectrum 
for m x = 100 GeV, T k a = 8 MeV and cross section values 
satisfying the constraint in the relic density: Vt x h 2 = 0.1143. 
The 2a observational upper limit on \/j,\ is 9 x 10 -5 . The plot 
shows the values of the logarithm of fi color-coded according 
to the scale on the right. 



The density of neutralinos changes with time (temper- 
ature through the variable x x ) according to the Boltz- 
mann equation. Its value at a given radiation temper- 
ature is related to the ratio Y given by Eq. [TT] (recall 
that p x = m x sY) by replacing xq with the corresponding 
x = m x /T and using the relation between the radiation 
temperature and that of the neutralino gas: 

Px( x x) = m x s(x = (x k dX x ) 1/2 )Y(x = (x k dX x ) 1/2 ) (16) 

Note that for the redshift range of interest T <470 eV 
and thus g*s, which enters in the value of the entropy 
s[x), is a constant equal to g*s( x o) — 3.915. 

Now, t kd (l + z kd ) 2 = t E (l + z E ) 2 = i (l + z,E) 1/2 , where 
ze is the redshift of matter and radiation equality and to 
is the age of the Universe. In the last equality we have 
used the relation between the scale factor and time for an 
Einstein-de Sitter Universe. Since (1 + Ze) = ^x,o/^7,0) 
we can finally write: 



P = 1-4/ 



Q x ^oh 2 



3/2 



^X,0Pcrit,0 



, ( Y{x={x kd x x f/ 2 ) 



t (<TV) S 
2 



dx x (17) 



Note that the function in parentheses inside the integral 
is equal to one if the comoving density of neutralinos is 
no longer evolving (has frozen) by z = 2 x 10 6 . This is 
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the case in the majority of the parameter space of the 
Yukawa interaction because by this redshift the velocity 
dispersion of dark matter particles is low enough for the 
Sommcrfcld enhancement to be saturated. As we men- 
tioned before, this means that the comoving density is 
essentially frozen afterwards. 

Fig. [3] shows the values of /i found by solving Eq. (fT"7|) 
for the same scan of the parameter space {m^/m x , a c ) 
as in previous figures. We took as fiducial values m x 
100 GeV, T kd = 8 MeV and il x h 2 = 0.1143, that is, the 
values of the normalization to the cross section are re- 
quired to give the observed dark matter abundance. The 
values of the logarithm of fi appear color-coded accord- 
ing to the scale in the right of the figure. The resonances 
inherited from S(x) and shown in Fig. Q] are also clear in 
this figure. Most of the region on the lower right corner 
is above, or very close to, the observational upper limit 
\fi\ < 9 x 10~ 5 . Points of the parameter space very near 
resonances also give values of /i that can be excluded. To 
show in more detail the behavior near a resonance, we 
present in Fig. [4] a plot of the values of fi as a function 
of the perpendicular distance, in parameter space, to the 
first resonance in the upper left of Fig. [3] Lines with dif- 
ferent colors show the magnitude of the change produced 
by different values of the kinetic decoupling temperature 
and of tt x h 2 , corresponding to the legends in the figure. 
The upper limit on /i according to the COBE/FIRAS 
experiment is also marked as a horizontal line. The fig- 
ure shows how close one can get to a resonance without 
violating the current constraint. 

This is a significant point, because values very close to 
a resonance are often invoked to provide a large boost 
and so to match observations of the positron excess in 
cosmic rays. 

According to Eq. (TTT1) . the value of fi decreases with 
increasing mass (/i cx \/m x ) because (av)s is nearly in- 
dependent of mass, and because, except very near reso- 
nances, S(x x ) is already saturated by z = 2 x 10 6 . In 
this way, our results can be easily scaled to the desired 
neutralino mass. Care is needed very close to the reso- 
nances since the simple scaling fi cx \jm x is no longer 
valid there. 

Finally, we show in Fig. [5] the analogous result for the 
values of fi in the case where the interaction responsible of 
the Sommerfeld enhancement is Coulomb-like (m^ — > 0) . 
In this case the enhancement depends only on a c . The 
plot shows, for the same fiducial case (m x — 100 GeV, 
T kd = 8 MeV and n x h 2 = 0.1143), that all cases with 
a c 6 x 10~ 2 are already excluded. 



[i, near first resonance 



COBE/FIRAS forbidden 




T kd =4 KT'GeV 

r td =8 l(T 3 GeV 

— r kd =i50io- r 

nh 2 =0.1075 

Oh 1 =0.1211 




-0.0015 -0.0010 -0.0005 0.0000 0.0005 0.0010 



"J 2 +log(m< t /m x ) 2 (pcrp . t0 



resonance) 



FIG. 4: The values of the fi distortion parameter near the 
first resonance in the upper left corner of Fig. [3] as a function 
of the perpendicular distance to the resonance. The lines 
with different colors correspond to different values of Tkd and 
£l x h 2 as marked in the legend. The solid black horizontal line 
shows the upper limit on n according to the COBE/FIRAS 
experiment. As in Fig. [3] m x — 100 GeV. 
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COBE/FIRAS allowed 
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SUMMARY AND CONCLUSIONS 



The prospects for dark matter detection have increased 
considerably in recent years. There is a continual ad- 
vance in the sensitivity of detectors on Earth that look 
for direct dark matter elastic scattering with nuclei, and 
in experiments that search indirectly for dark mater by 



FIG. 5: The same as Fig. [3] but for the case where the inter- 
action is strictly Coulomb-like. We have taken fiducial values 
of m x = 100 GeV, T kd = 8 MeV and fl x h 2 = 0.1143. Values 
above the black solid horizontal line are excluded by the cur- 
rent limits of /x according to the COBE/FIRAS experiment. 
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looking for non-gravitational signatures of the byprod- 
ucts of its hypothetical annihilation. Such technological 
improvements may lead, in the near future, to a definite 
proof of the existence of dark matter. 

Perhaps dark matter has produced non-gravitational 
signals that have already been detected. The recently re- 
ported excesses of positrons in cosmic rays by PAMELA 
and of electrons+positrons by FERMI/ ATIC/PPB- 
BETS can be explained by dark matter annihilation. Al- 
though other explanations with a different astrophysical 
origin are also possible, a solution based on dark matter 
is an attractive possibility. 

However, such solution seems to require an attractive 
force between the dark matter particles to enhance their 
annihilation through the Sommerfeld mechanism. For 
a Yukawa interaction via a single scalar, the magnitude 
of the enhancement depends on the coupling strength of 
the interaction a c , on the mass ratio of the force carrier 
to the dark matter particle m^/m x and on the relative 
velocities of the annihilation particles /3. 

Several recent papers have computed the boost to the 
annihilation due to a Sommerfeld enhancement (in the 
Galactic halo and/or in its subhalos) [ljl [2(| [Hf. One 
of the main aims of these works is to show that for cer- 
tain values of a c and m^/m x , this mechanism is able to 
produce large enough boosts to explain the cosmic ray 
anomalies. 

Only a handful of studies have addressed in detail the 
impact that such an enhancement has in the early Uni- 
verse. In [2l[, the authors found that if the cross sec- 
tion increases with decreasing relative velocity as 1/(3 
(which is valid in a certain regime for general Sommer- 
feld enhancement models), dark matter annihilation in 
the first halos would heat and ionize the IGM, violat- 
ing current constraints from the CMB. As was later 
pointed out by and [2(J, this problem is alleviated 
in more general models where the enhancement satu- 
rates at low velocities. In 22 1 and (23[, a constraint 
on the annihilation cross section was obtained by con- 
sidering limits on the energy deposition by annihilation 
at recombination. The constraint reported by [23j | is 
(tv)rec < 3.6 x lO-^cnrVVx/lTeVV/flEc, where 
Jb.ec is an average efficiency of energy injection into the 
IGM by annihilation at recombination. 

In the present paper, we have analyzed the impact of 
dark matter annihilation with Sommerfeld enhancement 
at higher redshift. 

In the first place, we have revisited the calculation 
of the dark matter particle abundance by solving the 
Boltzmann equation from freeze-out through the epoch 
of kinetic decoupling, including a full solution to the 
Schrodinger equation. Contrary to previous claims 
[20L H3] , we have found a significant suppression of the 
relic density, in agreement with a recent work by [3~j| . 
This suppression is particularly important near reso- 
nances, which are typically invoked to explain the cosmic 
ray anomalies. 

We found that to fit the observed dark matter abun- 



dance, the normalization of the cross section needs to be 
lowered by up to a factor of 10 compared to the case 
without enhancement (see Fig. [TJ). The result depends 
on the coupling strength and the proximity to a reso- 
nance. Exploring a broad range of dark matter particle 
masses and kinetic decoupling temperatures, we found 
a minor-to-medium impact of these parameters on our 
results; variations on T^d have the strongest impact. 

Secondly, we have calculated the amount of energy de- 
posited by dark matter annihilation into the radiation 
plasma in the redshift range 5.1 x 10 4 < z < 2 x 10 6 . 
Energy injection at this epoch would create a distortion 
in the CMB energy spectrum first pointed out by [28|. 
The COBE/FIRAS experiment has put constraints on 
this Bose-Einstein fi— type distortion. The upper limit 
on the chemical potential associated to the distortion is 
\fi\ < 9.0 x 10~ 5 at the 2a level [37[ . We have found 
that very near resonances, annihilation with Sommerfeld 
enhancement is already ruled out by this constraint on fi 
(see Fig [3]). Quantitatively, the "safe" range of proximity 
to a resonance depends on the particular model, higher 
values of dark matter mass and kinetic decoupling tem- 
perature allow a closer proximity to the resonance (see 
Fig. 2]). In the case where the force carrier is massless, 
the Yukawa interaction reduces to a Coulomb one. Our 
findings are much more stringent in this case with values 
ofa c > 6 x 10~ 2 already ruled out (see Fig. [5]). 

Improved upper limits on the fi parameter by a null dis- 
tortion detection in the CMB spectrum would rule out a 
larger region of the parameter space for the Yukawa inter- 
action. In [38], it was pointed out that an improvement 
close to two orders of magnitude has already been possi- 
ble for a number of years, and [4l| suggests that another 
order of magnitude is perhaps within reach. An upper 
limit on \fi\ of the order of 10~ 7 would certainly exclude 
large regions of the parameter space. It would exclude 
to a large extent near-resonance regions, which are the 
ones that produce the largest boosts to the annihilation. 
On the other hand, the detection of a distortion could 
possibly tell us something about the parameters of the 
Yukawa interaction and the nature of the force carrier <j>. 

In summary, our results indicate that for a given set of 
parameters (m x , T^d, m^, a c ), it is necessary to compute 
in detail the relic density to obtain the range of values for 
the normalization to the cross section that are compati- 
ble with current estimates of the dark matter abundance. 
Once this normalization is known, the energy input pro- 
ducing a fi— type distortion in the CMB spectrum can be 
computed to check whether the particular model violates 
current constraints. Only for allowed models can a boost 
factor for local dark matter annihilation be computed 
and advocated. 

Our findings show that the local boosts reported in the 
literature need to be renormalized to the proper value of 
(av)s implied by the observed relic density. This renor- 
malization can exceed a factor of 10 in extreme cases. 

A recent analysis of the PAMELA and FERMI data 
by [42J suggests that boost factors, over a standard value 
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BF o- = 0.0005 




FIG. 6: The values of the boost factor, relative to {av)s = 
3 x 10 _26 cm 3 s _1 for the parameter scan of the Yukawa interac- 
tion. These values are obtained according to the cross section 
values of Fig. [I] which are consistent with Q. x h 2 — 0.1143. A 
Maxwell-Boltzmann velocity distribution with a v = 5 x 10 -4 
for a neutralino with m x = lOOGeV and Tkd = 8 x 10 _3 GeV 
was used. The values are color-coded logarithmically accord- 
ing to the scale on the right. 

of ((Jv)s — 3 x 10 _26 cm 3 s _1 , of order 1000 or higher are 
required to explain both datasets simultaneously. We 
note that smaller boost factors (~ 100) can still fit the 
PAMELA data, but not the FERMI one, for a 100 GeV 
neutralino, which is the standard case we have considered 
here. These large boosts are also found in the analysis 
by [H[. It has been argued in the past that these boosts 
can be achieved only by invoking Sommerfeld-enhanced 
annihilation. For instance, [HJ and [33[ obtain maxi- 
mum boost factors ~ 1000 by assuming a local Maxwell- 
Boltzmann velocity distribution with a velocity disper- 
sion of 150 kms^ 1 (5 x 10~ 4 c), which roughly corresponds 
to the estimated local dark matter velocity dispersion. 
We find that these boosts are modified once the proper 
cross section is used. 

In Fig. [6] we show the boost factors, i.e., the mul- 
tiplicative factor to (av)s — 3 x 10 _26 cm 3 s _1 , that 
we find for a Maxwell-Boltzmann velocity distribution 
with <i v = 5 x 10~ 4 for the parameter space of the 
Yukawa interaction. The figure is for m x = lOOGeV, 



Tkd = 8 x 10 ~ 3 GeV and is consistent with a relic abun- 
dance Vt x h 2 = 0.1143. Keep in mind, however, that 
changes in the neutralino mass have almost no impact 
on the normalization of the cross section. Different ki- 
netic decoupling temperatures and variations on Q x h 2 
have also only a small effect (< 60% in combination). 
Thus, Fig. |5] is approximately correct for all the cases 
considered in this paper. The figure shows that even in 
extreme cases, for resonances in the upper right of the 
figure, the boost factors are < 500. In more favored re- 
gions of the parameter space a c > 10 -2 , m^/m x < 10~ 3 
(see for example (33|), the boost factors are < 100. 

This result suggests that additional assumptions are 
needed to account for the boost needed to explain the 
cosmic ray anomalies by dark matter annihilation alone. 
The inclusion of colder substructures to the overall 
smooth component with higher densities and lower ve- 
locity dispersions and thus, higher boost factors, could 
perhaps solve the issue. However, as found in recent 
high-resolution N-body simulations, the local dark mat- 
ter distribution is rather smooth 16]. Typical estimates 
on such an additional boost factor due to substructure 
in the galactic halo (albeit without Sommerfeld enhance- 
ment) are of the order of 1.4, unless a subhalo happens to 
be very close to Earth, in which case this boost could be 
larger than 10 These estimates do not include a pos- 
sible further amplification due to the scaling of the Som- 
merfeld enhancement with velocity dispersion. However, 
the constraints found by [23[ based on energy deposition 
at recombination suggest that the enhancement must be 
already close to saturation for a v = 150kms _1 , thus such 
a further amplification seems implausible. Therefore, is 
not clear that the inclusion of colder substructures can 
account for the additional boost. 

Finally, we mention that we have created a web 
application that solves the relevant equations de- 
scribed in this paper. It allows the user to compute 
individual Sommerfeld boosts, thermally averaged 
enhancements, cross section and /i values, among 
other quantities. The interested reader can find this at 
|http://www.m pa-garching.mpg.de / ~ vogelsma /sommerfeld / 
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